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ABSTRACT 

We have conducted a series of numerical experiments with the spherically symmetric, general relativistic, 
neutrino radiation hydrodynamics code A^;7e-B0LTZTRAN to examine the effects of several approximations 
used in multidimensional core-collapse supernova simulations. Our code permits us to examine the effects of 
these approximations quantitatively by removing, or substituting for, the pieces of supernova physics of inter- 
est. These approximations include: (1) using Newtonian versus general relativistic gravity, hydrodynamics, 
and transport; (2) using a reduced set of weak interactions, including the omission of non-isoenergetic neutrino 
scattering, versus the current state-of-the-art; and (3) omitting the velocity-dependent terms, or observer correc- 
tions, from the neutrino Boltzmann kinetic equation. We demonstrate that each of these changes has noticeable 
effects on the outcomes of our simulations. Of these, we find that the omission of observer corrections is par- 
ticularly detrimental to the potential for neutrino-driven explosions and exhibits a failure to conserve lepton 
number. Finally, we discuss the impact of these results on our understanding of current, and the requirements 
for future, multidimensional models. 

Subject headings: methods: numerical — neutrinos — radiative transfer — supernovae: general 



1. INTRODUCTION 



Colgate and White (1 19661) were the first to propose that 
core-collapse supernovae may be neutrino-driven and per- 
formed the first numerical simulations of such events, launch- 
ing more than four decades of research that continues to this 
day. A significant milestone occurred nearly two decades later 
with Wilson's discovery that delayed neutrino-driven explo- 
sions c ould be obtained. Based on his mode ls, Wilson con- 
cluded (IWilsonll 19851: iBethe & Wilsonlll985l) that the stalled 
supernova shock wave could be revived via neutrino absorp- 
tion on a time scale of several hundred milliseconds given the 
intense flux of neutrinos emerging from the proto-neutron star 
liberating the star's gravitational binding energy. Observa- 
tions of the neutrinos from S N1987A, the first such observa- 
tions of supernova neutrinos (iBionta et alJll987l : lHirata et alJ 
Il987h . provided support for the central role of neutrinos in 
the explosion mechanism. State-of-the-art simulations to- 
day continue to explore Wilson's neutrino-driven explosion 
mechanism in the context of two- and three-dimensiona l 
models (e.g., see Burrows et afl 120071; iMarek & Jankall2009l; 
iBruenn et alll2009^,Suwa et al1l2010l ). 

Neutrinos are weakly interacting particles whose cross sec- 
tions are energy dependent. Thus, unlike all other components 
in a supernova model, they are not well described as a fluid. 



|elen tz@utk.edu, mezzacappa a@oml.gov| 

' Department of Physics and Astronomy, University of Tennessee, 
Knoxville, TN 37996-1200, USA 

- Pliysics Division, Oak Ridge National Laboratory, P.O. Box 2008, Oak 
Ridge, TN 37831-6354, USA 

^ Joint Institute for Heavy Ion Research, Oak Ridge National Labora- 
tory, RO. Box 2008, Oak Ridge, TN 37831-6374, USA 

* Computer Science and Mathematics Division, Oak Ridge National 
Laboratory, RO.Box 2008, Oak Ridge, TN 37831-6164, USA 

^ National Center for Computational Sciences, Oak Ridge National Lab- 
oratory, RO. Box 2008, Oak Ridge, TN 37831-6164, USA 

' Department of Physics, University of Basel, Klingelbergstrasse 82, 
CH-4056 Basel, Switzerland 

' Department of Physics, Florida Atlantic University, 777 Glades Road, 
Boca Raton, FL 33431-0991, USA 



except in the deepest layers, and their transition in space to 
non-fluid-like behavior depends on their energy. Instead, the 
evolution of the neutrino radiation field, particularly in the 
semi-transparent regime, is far better characterized by classi- 
cal kinetics — specifically, the general relativistic Boltz mann 
kinetic equation (e.g., see lCardall & Mezzacappall2b03h . 
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Equation ([T]i describes the evolution of the neutrino distri- 
bution function f{t, xi,X2,X3, ni, E), which at time t 
and spatial location {xi,X2,xrj) supplies the distribution of 
neutrinos in direction cosines {111,112) and energy E'-i.e., 
the angular and spectral distribution of neutrinos. One such 
Boltzmann equation is solved for each flavor of neutrino — 
electron, muon, and tau {v^, v^, and Vt, respectively) — and 
for their antineutrino partners {v^, Vn, and v.,-)- The invariant 
colHsion term, C[/], in equation Q is written using emis- 
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sion, j, absorption, ™d scattering and pair kernels, R, 
foll owing the forms often used for neutrino transport (e.g., 
see ' Mezzacappa & Messeii[l99 9'), where / is the distribution 
function for the partner antineutrino and fi is the neutrino di- 
rection cosine. In equations ([T]i and (|2]l, / is a function of 
{fi, E), as well as position and time. The (/i', E') dependence 
of / and / inside the integrals illustrates the physical cou- 
pling of all energies and angles for each neutrino species and 
of neutrino and antineutrino partners. 

The first term on the left-hand side of equation ^ describes 
the time evolution of the local neutrino distribution owing to 
spatial transport through the volume of interest. The second, 
far more complex term on the left-hand side (LHS/Term 2) de- 
scribes the evolution of the local neutrino distribution in angle 
and energy as the result of (A) the coordinate system chosen, 
(B) special relativistic effects, and (C) general relativistic ef- 
fects. In what follows, we will refer to (B) as "observer cor- 
rections." 

Terms describing (A) depend on the choice of coordinate 
system. For example, in spherical-polar coordinates, the neu- 
trino direction cosine relative to the outwardly pointing radial 
vector changes as the neutrino propagates through a local vol- 
ume. This coordinate-system effect is included in LHS/Term 
2 and is present even in the absence of fluid motion or general 
relativity. For Cartesian coordinates, the neutrino direction 
cosines do not change as a result of the coordinate system 
choice alone and, consequently, such a term is absent. 

Terms describing (B) depend on the frame of reference 
chosen to measure the neutrino direction cosines and ener- 
gies. The comoving frame, with neutrino direction cosines 
and energies measured in an inertial frame of reference in- 
stantaneously comoving with the stellar core fluid with which 
the neutrinos interact, is often used. Neutrino-matter interac- 
tions are naturally expressed in this frame. Given this choice, 
the terms in LHS/Term 2 present a significant numerical chal- 
lenge. Finding discrete representations that guarantee conser- 
vation of lepton number and energy is one of the most difficult 
aspects of modeling neutrino transport in stellar cores. This 
has been achieved for general relativistic, spherically symmet- 
ric flows (Liebendorferet al. 2004), providing the conceptual 
and implementation groundwork for achieving the same in ax- 
isymmetric (2D) and non-symmetric (3D) flows. Further the- 
oretical foundations have been laid (ICardall & Mezzacappal 
I2003t ICardaU et al]l2005h : steps toward the development of 
a 2D Boltzmann solver have been taken (Ott et al. 2008); and 
the challenge now is to fully implement lepton energy and 
number conserving discretizations in 2D and 3D models. 

For another choice of reference frame — measuring neutrino 
angles and energies relative to the inertial, "lab" frame of 
a distant observer — the terms encapsulating the special rel- 
ativistic effects in LHS/Term 2 are absent, simplifying the 
left-hand side of the Boltzmann equation. In such a frame of 
reference, the neutrino direction cosines and energies do not 
change from observer to observer in the frame. However, this 
simplification comes at a price because the neutrino-matter 
interactions are naturally described in the comoving frame. 
In the lab frame, a Lorentz transformation is required in or- 
der to express the comoving-frame neutrino-matter interac- 
tions in terms of the lab-frame direction cosines and energies, 
which introduces non-trivial velocity dependencies into the 
lab-frame collision term. 

One approach to the complexity of the lab-frame colli- 
sion term is the "mixed frame" approach, which uses the 



lab-frame 4-momenta and an 0{v/c) Taylor-series expan- 
sion i n energy of the comovin g- frame emissivities an d opac- 
ities ( iMihalas & Kleiiil 11982b . iHubenv & Burrowsl (l2007h 
have proposed to use the mixed-frame approach for core- 
collapse simulations with extensions for non-isotropic and 
non-isoenergetic scattering. The mixed-frame approach may 
be difficult to extend to arbitrarily relativistic flows, and has 
not yet been used in the context of a full-physics core-collapse 
supernova simulation. 

In a general relativistic setting, such as core-collapse su- 
pernovae, we must contend with (B) and/or (C) regardless of 
the frame of reference chosen to describe the neutrino direc- 
tion cosines and energies. Even for static general relativistic 
environments, angular aberration, gravitational red shift, and 
other effects occur, and the resulting terms in (C) are always 
present. 

Regardless of approach, comoving- or lab-frame, it is prob- 
lematic to adapt the simplicity of both approaches, simultane- 
ously simplifying the left- and right-h and sides of the Boltz- 
mann equation, as has been done in iBurrows et al.l (120061 
I2007h . lOtt et al.1 (120081) . and other models using die Vul- 
can/2D code, although one can view the implementation in 
these works as steps toward a more complete description. 
They deploy a lab-frame approach for terms describing angu- 
lar aberration and energy shift on the left-hand side (or assume 
such terms are unimportant in a comoving-frame approach), 
while simultaneously deploying a comoving-frame approach 
for the collision term describing the neutrino-matter interac- 
tions on the right-hand side. This is not a mixed-frame ap- 
proach in the sense described above. It is an approach not 
based in any reference frame, and it is physical only for static 
cases in which there is no distinction between lab and comov- 
ing frames. One of the goals of this study is to investigate 
the importance of the terms in LHS/Term 2 in a comoving- 
frame approach, and whether they can be ignored while using 
a comoving-frame approach for the collision term. 

Modeling general relativistic Boltzmann kinetics is also 
challenging because of the complexity of the collision term 
on the right-hand side of the Boltzmann equation, even in 
a comoving-frame formulation. Looking at equation (O, 
we see that the collision term describes the full, direct cou- 
pling of all neutrino angles and energies for each neutrino 
species, owing to neutrino isoenergetic (IS) scattering on nu- 
clei and non-isoenergetic (NIS) scattering on electrons and 
nucleons. The pair creation and annihilation processes (PR) 
such as electron-positron annihilation and nucleon-nucleon 
bremsstrahlung also couple the angles and energies of the neu- 
trino and antineutrino species of each flavor together The 
coupling of all neutrino angles and energies through the rele- 
vant set of weak interactions dominates the computation asso- 
ciated with the solution of the neutrino Boltzmann equations. 
It has been argued (Burrows et al. 2006, 2007; Nordhaus et al] 
2010) these couplings are subdominant and can be ignored, 
greatly simplifying the neutrino Boltzmann equations and sig- 
nificantly reducing the computational cost associated with 
their solution. A second goal of this study is to investigate 
whether or not such approximations to the collision term are 
realistic. 

The complete general relativistic Boltzmann equation 
was solved in spherically symmetric models of core- 
collapse supernova e by the Oak Ridge-Basel collaboration 
dLiebendorfer et aLllTO Ol, 2004) and by Sumiyoshi and col- 
laborators ( ISumivoshi et al.ll2005() . Achieving this in three- 
dimensional models of core-collapse supernovae presents a 
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major challenge, one that will likely require sustained exas- 
cale resources to meet. 

The overarching goal of this study is to use gen- 
eral relativistic, spherically symmetric Boltzmann simula- 
tions to guide and, more importantly, set minimum re- 
quirements for accurate 2D and 3D simulations. We 
use the Oak Ridge-Basel code Ag/Ze-BOLTZTRAN in these 
studies to compare general relativistic-full weak interac- 
tion physics (GR-FullOp), Newtonian-full weak interaction 
physics (N-FullOp), Newtonian-reduced weak interaction 
physics (N-ReducOp), and Newtonian-reduced weak inter- 
action physics-no observer correction (N-ReducOp-NOC) 
models. These models will demonstrate the importance of 
general relativity, a complete weak interaction set and treat- 
ment, and the terms in LHS/Term 2 to stellar core collapse 
and the post-core-bounce evolution. Current multidimen- 
sional models suggest that spherical symmetry is a reason- 
able approxima tion for t he first 100-150 ms after bounc e 
jMarek & Janka 2009; Br uenn et al.ll2009HSuwa et al.ll20Tol) . 
Thus, the simulations presented here are relevant for dis- 
cussing the initial conditions present for all multidimensional 
phenomena that might ensue; e.g., neutrino-driven convection 
and the standing accretion shock instability (SASI). 

2. DISABLING OBSERVER CORRECTIONS IN A 
LAGRANGIAN FORMULATION 

Disabling general relativity in a simulation, instead running 
a Newtonian simulation, is straightforward and requires no 
special considerations to define or interpret. The same holds 
true for limiting the weak interaction channels included in the 
collision term on the right hand side of the Boltzmann equa- 
tion. Newtonian or general relativistic simulations can be per- 
formed with more, or less, weak interaction physics. How- 
ever, disabling the observer corrections in a model requires 
some de finition and care. 

Using iMezzacappa & Matzne^ (11989 ). equation (VI. 11), 
we begin by expressing the neutrino Boltzmann equation in 
flat spacetime, Eulerian spherical-polar spacetime coordinates 
with zero shift vector, comoving-frame 4-momenta, and in 
nonconservative form: 
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In equation (O, mo and Eq are the neutrino direction cosine 
and energy as measured in a comoving frame of reference, 
and e and o are the invariant emissivity and opacity. We use 
c = 1 throughout this section, and have written the Eulerian 
time coordinate as i. Multiplying by (l+fiov) and rearranging 
we have 
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We can express the observer correction terms in conservative 
form using 
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Substituting these expressions into equation (|5]l, we get the 
0{v) Boltzmann equation. 
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Note, the last term on the left hand side of equation (|6]), 
— {f/p)Dp/Dt, is part of the observer corrections. 

Our Eulerian starting point with comoving-frame neutrino 
4-momenta (equation lU provides the formulation in which 
we can most readily discuss what it means to have Boltzmann 
neutrino transport without observer corrections. We define 
the evolution of the neutrino distributions in this case to be 
governed by 
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that is, by equation dUl with the velocity-dependent, neutrino 

angle- and energy-shift terms ign ored. 

From equations (15)-(22) of iLiebendorfer et al] (120041) . 
Ag;7e-B0LTZTRAN evolves the following purely Lagrangian 
(comoving-frame spacetime coordinates and neutrino 4- 
momenta) equations in flat spacetime. 
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where F — f /pis the specific neutrino distribution function. 
Note, here the Lagrangian partial derivatives with respect to t 
and m are at constant m and t, respectively. We can express 
the partial derivatives with respect to the Lagrangian space- 
time coordinates, {t, m), in terms of the Eulerian spacetime 
coordinates, {t,r), using 



and 



Using equation (45) of iCardall & Mezzacappal (I2003h . we 
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Substituting these transformations into the first three terms in 
equation ds} gives 
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Writing the right hand side of equation (|9]l in terms of / and 
expanding, we are left with 
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Moving the term containing the Lagrangian time-derivative, 
D/Dt, of the density to the LHS and restating it using the 
Lagrangian time-derivative, d/dt, we have to 0{v), 
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Therefore, with equation (|7J as a guide, a "no-observer- 
correction" run in our Lagrangian formulation would corre- 
spond to a solution of the following equation: 
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Equation ( fT2] i is used in our 'no-observer-corrections' model, 
N-ReducOp-NOC. Equation ( fT2b is clearly not manifestly 
conservative for neutrino number when integrated over mass 
and when the density evolves, and therefore its Eulerian 
equivalent, equation (|7]), is also not number conservative 
when integrated over volume. On the other hand, equation 
^ is manifestly number conservative when integrated over 
volume (after dropping the novdf/dt term). The culprit in 
equation ^ is the vdf /dr term. When expressed in volume- 
conservative form, this term contains a velocity divergence, or 
equivalently a logarithmic time derivative of the density, that 
would normally be cancelled by the logarithmic time deriva- 
tive of the density in equation But when the observer 
corrections are dropped, the last term on the LHS of equation 
(|6| is dropped, and this cancellation no longer occurs and we 
are left with the same term that appears as the last term on the 
LHS of equation (fTSl) . which breaks number conservation. By 
expressing our observer corrections in equation (|6]l in conser- 
vative form, we made explicit this logarithmic time derivative 
of density contained within them. 

3. NUMERICAL METHODS AND INPUTS 

All models in this paper are computed using the paral- 
lel version of the general relativistic, spherically symmetric, 
neutrino radiation hydrodynamics code A^//e-BOLTZTRAN 
(Lieb endorfer et al.. ,2004.) with extensions that we describe 
here. 

3.L Agi7e-B0LTZTRAN 

Ag//e-BOLTZTRAN is a combination o f the general rela- 
tivisti c (GR) hydrodynamics code Agile dLiebendorfer et all 
120021) and the neutrino transport code BOLTZTRAN 
jMezzacappa & Bruenn 1 993"btlMez zacappa & Messer* ! 19991: 
iLiebendorfer et al. 2004). Agile solves the complete GR 
spacetime and hydrodynamics equations implicitly in spher- 
ical symmetry on a dynamic, moving grid. The moving 
grid allows adequate resolution of the shock using only 
0(100) radial zones. Recent enhancements include the use 
of a TVD (total variation diminishing) hydrodynamics solver 
(Liebendorfer et al.1 12005). which improves the accuracy of 
advection, and the use of dm as the gri d coordinate rather than 
the enclosed mass dFischer et al.l201d §2.1), which improves 
numerical accuracy when mass zones are small and density 
gradients are large. In Newtonian mode the gravitational mass 
is set equal to the baryonic mass (omitting the non-rest-mass 
energy contributions) and the relativistic parameters are set to 
their non-relativistic values: a = l,r = 1. BOLTZTRAN 
(iMezzacap pa & Bruennlll993btlMezzacappa & Messedll999l: 
Liebendorf er et al.l 120041) solves the GR extension of the 
spectral neutrino Boltzmann equation (Eq. [8]l with a Gauss- 
Legendre (Sn) quadrature. Here we use an 8-point angu- 
lar quadrature and 20 logarithmically-spaced energy groups 
wit h group centers from 3 to 300 MeV. Previous stud- 
ies (IMezzacappa & Bruen"nlll993at ILiebendorfer et alj|2004h 
with (Ag;7e-)B0LTZTRAN have shown that 20-group energy 
resolution is adequate in removing artifacts seen at lower (12- 
group) energy resolution, and their 12- and 20-group runs ex- 
hibited no differences in outcomes. Moreover, 20-group en- 
ergy resolution matches, or exceeds, the resolution used for 
supernova models computed with the multidimensional codes 
we discuss in §??. The discretization scheme is designed 
to simultaneously conserve lepton number and energy as de- 
scribed in .Liebendorfer et al.. (.2004) . Since we do not include 
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Table 1 

Neutrino Opacity Summary Table 
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any physics to distinguish between muon- and tau-flavored 
leptons, we use the combined species v^r — {v^,Vr} and 

For all models we use the nuclear, ele ctron, and pho- 
ton equations of state (EoS) of Lattimer & Swestvl (Il991h 
with the bulk incompressibility of nuclear matter = 
220 MeV0 This match es the current experimental value of 
= 240 ± 20 MeV (iShlomo et al.1 12006) better than the 
value of 180 MeV more commonly used with LS EoS in 
the past, though the value of k,s in LS EoS has been shown 
to be of little consequence during the ea rly phases of core- 
collapse supernova evolution shown here (iSwestv et alJI 19941 
[Thompson et al.ll2003t Lentz et al.ll2010h . Matter outside the 
"iron" core is treated as an ideal gas of ^^Si that "flashes" in- 
stantaneously to nuclear statistical equilibrium when the tem- 
perature exceeds 0.47 MeV. 

The stellar progenitor used for all mode ls reported here is 
the \5 -M(^ solar-metalicity progenitor of IWooslev & Hege3 
(l2007h . We have mapped the inner 1.8 Mq of the progenitor 
onto 108 mass shells of the adaptive radial grid. 

3.2. Neutrino Opacities 

The base, or full, opacity set (FullOp) inc ludes emis- 
sion, absorption, and scattering on free nucleons (iReddv et al.l 
Il998h : isoenergetic scattering on a-particles and heavy 
nuclei ('Br uennI 1198 5): scattering of neutrinos on elec- 
tron s (NES) and positrons (NPS) (Schinder & Shapiro 
Il982h : production of neutr i no p airs from e+e^ anni- 
hilation (ISchinder & Shapir^ '1982) and nucleon-nucleon 
bremsstrahlung jHan nestad & Raffelt 1998); and electron 
capture (EC) on nuclei using the LMSH EC table of 
Langanke et al. (2003), which u tilizes the EC rates of 
Langanke & Martmez-Pinedd ^2000). The full angle and en- 
ergy exchange for scattering between the neutrinos and elec- 
trons, positrons, and nucleons is included, while scattering on 
nuclei is isoenergetic (IS). Bremsstrahlung and e^e+ annihi- 
lation are the only sources of Vf^^r and P^t-. 

For our reduced opacity set (ReducOp) we replace the 
LMSH EC table for electron capture on nuclei with an in- 
dependent particle approximation (IPA) (Fuller 1982) us- 
ing the implementation described in B ruenn (1985), which 
cuts off when the mean neutron number of the heavy nu- 
clei N > 40. We also drop all scatterings (NIS) that couple 
neutrino-energy groups. The primary contribution of electron 

^ We use the latest version of the ILattimer & Swestvl 119911) EoS, 
version 2.7, which is available for download from its authors at 
|http : / / www .astro ■ suny sb ■ edu/dswesty/lseos ■ html| 



and positron sca ttering opacities is through ne utrino-energy 
down-scattering dMezzacappa & Bruennll 1993d) . not through 
their contribution to the total scattering opacity; therefore, we 
omit them completely from the ReducOp opacity set. We also 
replace the NIS nucleon scattering opacities of iReddy et al. 
( 1^998) with the more approximate IS equivalents from lBruenn 
( Il985i) . For consistency, we al so replace the neut rino emis- 
sion and absorption opacities of lReddv et al.l (11998b with their 
iBruennI (|1985) equivalents. Ion-ion correlations and weak 
magnetism are omitted from both opacity sets. The two opac- 
ity sets are summarized in Table [T] 

3.3. Observer Corrections 

As noted in Sj2] the Lagrangian formulation in Agile- 
BOLTZTRAN and the use of the specific neutrino distri- 
bution function, F — f / p, which is needed to properly 
account for number and energy conservation (see discus- 
sion on the necessity of usin g F for Lagrangian models in 
ICardall & Mezzacappall2003l §IV.B), require care in the def- 
inition of a no-observer-corrections model. Moreover, time 
derivatives at fixed Lagrangian mass coordinates on a moving 
grid must be handled with care (see Liebendorfer et al. 2004 
§3.2). Therefore, for model N-ReducOp-NOC, we implement 
the "compression" term in the no-observer-corrections trans- 
port equation (fT2] i by re-expressing the time derivative of den- 
sity as a spatial divergence, using the continuity equation, 

Fdp _^ F d {rH) ^^^^ 
p dt dr 
4. RESULTS 

We present results from four spherically symmetric, core- 
collapse supernova models of decreasing physical fidelity. 
The most physically complete model (GR-FullOp, black lines 
in plots) utilizes the more modern and complete FuUOp opac- 
ities and the full general relativistic treatment of gravity, hy- 
drodynamics, and transport as described in Lieben dSrfer et al.l 
(2004 ). The first approximate model (N-FullOp, red lines) 
replaces the full general relativity of the GR-FullOp model 
with Newtonian gravity, 0{v/c) hydrodynamics, and 0{v/c) 
transport. The second approximate model (N-ReducOp, green 
lines) further replaces the more complete FullOp neutrino 
opacity set with the ReducOp opacities (see Table [T]for a full 
comparison), while retaining the Newtonian gravity, 0{v/c) 
hydrodynamics, and 0{v/c) transport. This approximation 
includes the important effect of removing neutrino weak in- 
teractions that down-scatter the neutrino energy. The final ap- 
proximate model (N-ReducOp-NOC, blue fines) retains the 
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Figure 1. Properties of our models at core bounce, where bounce is defined as the maximum compression of the central density during the launch of the bounce 
shock. The models are: general relativistic gravity, hydrodynamics and transport with full opacities (GR-FullOp, plotted in black); Newtonian gravity with 
full opacities and 0{v/c) hydrodynamics and transport (N-FullOp, plotted in red); Newtonian gravity with reduced opacities and 0{v/c) hydrodynamics and 
transport (N-ReducOp, plotted in green); and Newtonian gravity with 0{v/c) hydrodynamics and reduced opacities, and 0(1) transport (N-ReducOp-NOC, 
plotted in blue). The panels are: radial velocity (upper left), density (upper center), entropy (upper right), temperature (kT, lower left), net electron (or proton) 
fraction (Ye, lower center, solid lines), net lepton fraction (Y^ = Ye + (n^^ — ni7e)/™baryons. lower center, dashed lines), and pressure (lower right). All 
quantities are plotted relative to enclosed rest-mass in Mq . 



Table 2 

Model Approximations and Properties 



Property 


GR-FullOp 


N-FullOp 


N-ReducOp 


N-ReducOp-NOC 


Gravity and hydrodynamics 


GR 


Newtonian 


Newtonian 


Newtonian 


Neutrino opacities (see Table[T) 


Full 


Full 


Reduced 


Reduced 


Observer corrections 


GR 


Oiv/c) 


Oiv/c) 


None 


Homologous core, Mgii (Mq) 


0.429 


0.492 


0.717 


0.427 


Central-core density at bounce, pc (x lO^** g cm^^) 


4.714 


4.264 


3.336 


3.157 


Central-core electron fraction (Ve ) at bounce 


0.2448 


0.2407 


0.3046 


0.1855 


Central-core lepton fraction (Yl) aX bounce 


0.2804 


0.2782 


0.3696 


0.2007 


Peak shock radius (km) 


162 


190 


171 


142 


Peak -Luminosity (Bethe s~^) 


406 


450 


448 


160 



Newtonian gravity and 0{v/c) hydrodynamics of the previ- 
ous model, but drops the observer corrections completely, re- 
ducing the transport to [N.B. The and 0{v/c) hy- 
drodynamics equations are identical.] The general and core- 
bounce properties of all models are summarized in Table |2] 
and plotted in Figure [T] 

4. 1 . GR versus Newtonian Gravity 

The effects of general relativity on the core dynamics are 
seen in the comparison of the first two models (GR-FullOp 
and N-FullOp). The deeper gravitational well of the GR 
model results in a more compact homologous core at bounce 
(0.429 Mq versus 0.492 A/©) with a higher central density 
(4.71 X 10^'^g cm^'^ versus 4.26 x 10^'' g cm^'^) and higher 
temperatures throughout the unshocked core (Figure [U. The 



electron (Ye) and lepton (Yl) fractions are essentially un- 
changed modulo the shift in shock position, as are the velocity 
and entropy, while the pressure differences follow the density 
differences. As the shock moves out, the shock radius for 
both models (Figure |2]) remains close for the first 40 ms after 
bounce and then diverges. The GR-FullOp model has max- 
imum shock extent that is 30 km (20%) smaller than the N- 
FullOp model, and by 150 ms after bounce the shock radius is 
40 km (30%) smaller. Several quantities reflect the long-term 
effect of the more compact, and therefore hotter, proto-NS 
in the GR-FullOp model, including the higher luminosities 
for all neutrino species (Figure O and the higher RMS ener- 
gies ((i?,y)RMs) of neutrinos of all flavors after the break-out 
burst (Figure T hese differences are in acc o rd with those 
already reported by iLiebendorfer et al] ( l260l]) . iBruenn et alj 
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Figure 2. Shock trajectories in km, versus time after bounce, for all models. 
The colors have the same meaning as in Figure[T] Shock position is computed 
by bisecting the pair of mass shells with the largest negative radial velocity 
gradient —dvr/dr. 
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Figure 3. Comoving-frame neutrino luminosities measured at 400 km for 
all models. Colors are as in Figure [T] Electron neutrino, Ue, luminosities 
are represented by solid lines, Pe -luminosities by dotted lines, and j/^t- 
luminosities by dashed lines. -luminosities are indistinguishable from 
-luminosities, and omitted from this figure. The luminosities are in 



Bethe s , where 1 Bethe = 10 ergs. The lower panel provides a detailed 
view of the luminosities below 40 Bethe s~^ during the first 100 ms after 
bounce. 



(|2001'), and iBuras et al.1 (200^ using different progenitors, 
different opacity sets (similar to ReducOp, though including 
NES), different energy and angle resolutions, and for the latter 
two cases, different codes. Our GR/Newtonian comparison is 
included here for completeness and to facilitate relative com- 
parisons across all four models. 

4.2. Reduced Neutrino Opacities 
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Figure 4. Comoving-frame neutrino RIMS energies, (i?jy)RMS = 
(/ dfidEE^F/ J dfidE E'^Py/'^, measured at 400 km for all models. 
RMS energy is computed over number density^ not number flux. Colors are 
as in Figure [T] Line styles are as in Figure [5] The lower panel provides 
a detailed view of (i?,y)H,MS for values less than 20 IVIeV over the perious 
±40 ms. 



The changes induced as we go from the FullOp opaci- 
ties (model N-FullOp) to the ReducOp opacities (model N- 
ReducOp) in the Newtonian-gravity, ©(w/c) -hydrodynamics, 
and C'(ti/c)-transport limit are more dramatic than those seen 
for t he transition from models GR-FullOp to N-FullOp in 
i )4.1l The shock position at bounce changes from 0.492 Mq 
for N-FullOp to 0.717 Af© for N-ReducOp (Figure [B, with 
the entropy peak (upper right) making the same shift. The in- 
crease in the initial shock mass, Mgh, is correlated with the 
corresponding increase in core lepton fraction, from Yl = 
0.28 to 0.37 (Msb oc Yl). The larger M^h for N-ReducOp, 
relative to the other models, results in a correspondingly larger 
region of high pressure, temperature, and density at bounce. 
The vigorous post-bounce shock of model N- ReducOp results 
in a strong "ringing" of the shock (Figure|2|. iThompson et alj 
(2003) reported a similar ringing for their "no NES" model. 

The -luminosity of the N-ReducOp model reaches the 
same peak value as in the N-FullOp model, 450 Bethe s^^, 
but the breakout burst is much shorter in duration and 
represents a smaller total emission of ly^- The shock 
starts out at a larger mass coordinate and passes through 
less total mass before becoming a steady accretion shock. 
Like Thompson et al. ( 2003), we see oscillations of the i^- 
luminosities and (i?j/)RMS (Figures [3] & S]) just after bounce 
induced by shock oscillations passing through the neutri- 
nospheres. 

The differences between the N-FullOp and N-ReducOp 
models can be understood by considering three opacity 
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changes imposed simultaneously: (A) the inclusion of 
NES/NPS; (B) the use of th e LMSH EC table; and (C) the 
use of the Reddv et alj ( Il998h nucleon opacities. 

(A) The effects of omitti ng the NES opacity alone durin g 
collapse were explored by Mezzacappa & BruennI (Il993cl) . 
who showed that energy down-scattering by NES allowed the 
energy down-scattered neutrinos to escape more easily be- 
cause of the lower absorption and scattering cross-sections at 
lower energies, and reduced the core Ye by 15% and neutrino 
fraction, Y,y, by 30% in their model with NES relative to one 
without NES. The higher number of trapped neutrinos, with- 
out NES, is reflected in the higher core Ye, Y ^, and Yi, for our 
model N-ReducOp. The "no-NES" model of [Thompson et al] 
(I2OOI (§7.4, Fi gures 20 & 21) also shows large differences in 
(£'y^^^}FiMS> with a bounce "spike" reaching 32 MeV, and a 
{E^^^)-[{MS increase at 150 ms post-bounce of 4 MeV (20%) 
relative to a model with NES. This compares to a 50 MeV 
"spike" and 7 MeV (40%) increase at 150 ms after bounce 
in (£'i/^^)rms for our N-ReducOp model relative to our N- 
FullOp model. 

(B) We have removed the LMSH EC table from ReducOp 
opacities in favor of the simpler IPA prescription, as not all 
modern supernova simulations use EC rates like the LMSH 
EC table, which reflect the ense mble of nuclei an d their ex- 
cited states in the collapsing core. iHix et all (12003b found that 
the enhanced EC arising from the removal of the artificial cut- 
off in IPA for heavier nuclei that occur at higher densities dur- 
ing collapse decreased the central-core Ye and Msh at bounce 
by 10% and 20%, respectively, relative to the IPA implemen- 
tation. Conversely, IPA overestimates EC where it is active 
and leads to excess deleptonization and stronger collapse in 
the outer regions of the Fe-core where the > 40 cut-off cri- 
terion is not triggered. One such region can be seen at bounce 
outside the homologous core near 0.9 Mq (Figure [T]i, where 
lower Ye and higher density exist in the N-ReducOp model 
relative to the GR-FullOp and N-FullOp models. 

(C) We have also replaced the i/-nucleon emission, absorp- 
tion, and scattering opacities (Reddy et al. 1998) in FullOp 
with the corresponding opacities of iBruenni (11985 ) to elimi- 
nate neutrino-ene rgy down-sc attering on nucleons. We have 
previously found dLentz et al.i l2010) that the inclusion of the 
enhanced nucleon opacities results in an enhancement of the 
luminosities, but not RMS energies, and lifts the post-bounce 
shock outward by 10 km by 100 ms post-bounce through ab- 
sorption of the excess luminosity. These findings are consis- 
tent with the results of R ampp et al.l 020021) on the enhanced 
neutrino-nucleon opacities. 

4.3. No Observer Corrections 

For the final comparison we change the treatment of the 
observer corrections in the transport equation. In model N- 
ReducOp-NOC we have removed the velocity-dependent ob- 
server corrections from the Boltzmann transport equation in 
the appropriate, Lagrangian approach as described by fj2] 
equation (fT2] i. but retain the Newtonian gravity, 0{v/c)- 
hydrodynamics, and reduced opacities of the N-ReducOp 
model. 

Dropping the observer corrections in model N-ReducOp- 
NOC results in a dramatic change in the properties of the core 
at bounce, as can be seen in Figure[T]and Table|2] The homol- 
ogous core mass drops from 0.717 Mq for the N-ReducOp 
model to 0.427 Mq for the N-ReducOp-NOC model; the lat- 
ter of which is virtually indistinguishable from the homolo- 
gous core mass, A/gh, for the most physically complete model 



(0.429 Mq for GR-FullOp). This coincidental alignment of 
the bounce shock positions for the most and least physically 
complete models should be contrasted with the lower electron 
and lepton fractions (Ye, Yl) and density (see Table |2]i in the 
homologous core for the N-ReducOp-NOC model relative to 
the GR-FullOp model, as well as the larger inflow velocities 
and densities and lower Ye and Yj, outside the shock, which 
imply, among other things, an increased ram pressure against 
which the shock must propagate. 

The shock (Figure |2]i in the N-ReducOp-NOC model starts 
more vigorously than in the FullOp models, but does not show 
the large overshoot of the other reduced opacity model, N- 
ReducOp. All of the shock trajectories cross near 35 ms after 
bounce, with the N-ReducOp-NOC model having the deepest 
shock throughout the rest of the run. This is in stark contrast 
to the other Newtonian models, which have larger shock radii 
relative to the GR-FullOp model. 

The neutrino luminosities (Figure [3]i of the N-ReducOp- 
NOC model are also substantially lower than for any other 
model. The i^e-luminosity from shock-breakout peaks at 
160 Bethe s"^ for the N-ReducOp-NOC model relative to 
the 400^50 Bethe s^^ for the models with observer correc- 
tions. The Pg-luminosities of all models approaches the Ve- 
luminosities at around 80 ms after bounce, and track together 
thereafter By 150 ms after bounce, the i/gPe-luminosities for 
the N-ReducOp-NOC model are approximately 32 Bethe s^^, 
while the other two Newtonian models have Ve Pe-luminosities 
of ^ 55 Bethe s~^ and the GR-FullOp model has Vei^e- 
luminosities of ^ 60 Bethe s^^. 

The RMS neutrino energies (Figure |4| for h'e (solid lines) 
and De (dotted lines) in the N-ReducOp-NOC model fol- 
low those in the N-ReducOp model closely, with (£',y^)RMS 
slightly higher before bounce and at later times, except in 
the immediate post-bounce period when it oscillates in the 
N-ReducOp model with the shock. In contrast, the post- 
bounce {E^^^)-[{MS for model N-ReducOp-NOC is substan- 
tially lower relative to model N-ReducOp, and the "spike" 
after bounce is gone. After breakout both ReducOp models 
have large (-E,y^^}RMS relative to the FullOp-models. 

As we traverse our four models, it is clear that all of the neu- 
trino luminosities are significantly affected. The general trend 
is for the luminosities to decrease considerably as we go from 
GR-FullOp to N-FullOp to N-ReducOp to N-ReducOp-NOC. 
The largest variations among the models are exhibited by the 
electron-flavor neutrinos, with luminosity variations as large 
as 35 Bethe s"^ at 150 ms after bounce. However, the vari- 
ations in the {E^Jj^us and (-EpJrms are not as dramatic as 
we traverse the four models, and not monotonically decreas- 
ing as the model sophistication decreases. Variations in the 
)rms^ like their luminosity counterparts, remain signifi- 
cant, although not monotonically decreasing with model. The 
post-bounce RMS energies, (i?iy^^)RMSi vary by > 10 MeV 
up to 150 ms after bounce. 

At 100 ms after bounce (Figure |5]) the shock encloses 
1.45 A/0 for all four models. The most significant differ- 
ences seen in the dense core are in Ye and (lower center), 
where the N-ReducOp-NOC model has generally the lowest 
Ye and Yl and lacks the peak the other models exhibit just 
outside their bounce shock positions, Msh- Figure |6] shows 
the same data, but focuses on the outer, "hot-mantle" region 
between the shock and proto-NS with the most evident dif- 
ferences being those related to the shock radius (smallest for 
N-ReducOp-NOC). Though shifted in radius, the hot-mantle 
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Figure 5. Same as in Figure □ but at 100 ms after core bounce. 
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Figure 6. Same as in Figure|5] but as a function of radial coordinate, r, in km. Net lepton number has been omitted, as Yj, for all but the inner few km. 
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Figure 7. Lepton number on the grid plus neutrino flux through the outer 
boundary, A'^L.cons, indicating the quality of numerical lepton conservation 
for all models. Model N-ReducOp-NOC loses 34.7% of the original AfL,cons 
by 200 ms after bounce, while the other three models conserve A'^L.cons to 
within 0.1%. 

region is similar among all models. 

The N-ReducOp-NOC model shows (Figure |2ll a drop in 
total conserved lepton number (iVL.cons, the lepton number 
on the computational grid plus the time-integrated number 
flux of neutrinos at the outer boundary) starting just before 
bounce and continuing throughout the rest of the run. This 
is not seen in the other models, which maintain lepton con- 
servation. The root of the non-conservation can be seen in 
the integration of equation (fTSI ) for neutrino number over the 
entire grid. It results from the {F / p){d p / dt) "compression" 
term and is strongest during the epoch of high and rapid 
density changes surrounding core bounce. 

We illustrate this loss through (Figure[8] solid lines) and 
Yl (dashed lines) profiles for the GR-FullOp model (upper 
panel) and the N-ReducOp-NOC model (lower panel) as a 
series of temporal snapshots near bounce. Until the shock 
reaches 0.5 M0 in the GR-FullOp model (upper panel), the 
core is still opaque, and the neutrinos are trapped. Therefore, 
total lepton number is conserved locally, and Y^ is steady in- 
side 0.5 Mq during shock breakout as depicted in the upper 
panel of Figure[8] When the shock reaches 0.5 Mq, it begins 
to break through the neutrinospheres, and the neutrinos can 
escape, causing the local Yl and Yf. to drop behind the shock. 
The escaping neutrinos contribute to Yl in front of the shock 
as a visible pulse, a small portion of which are absorbed by the 
cold, infalling matter ahead of the shock, forming a transient 
radiative precursor in Y^. For the N-ReducOp-NOC model 
(lower panel) there is no such corresponding epoch of local 
lepton conservation as the shock forms in the opaque (neu- 
trino trapped) inner core before emerging through the neutri- 
nospheres. The net effect of the compression term in equa- 
tion ( fT2] i is one of destroying neutrinos, which then results 
in a net decrease in Y^ via the interactions e^p o v^n and 
e^n DeP- The loss of neutrinos to the compression term 
reduces the neutrino pulse ahead of the shock and lowers the 
z/e-luminosity in the breakout burst (Figure[3]l. 

The fundamentally different behavior of the N-ReducOp- 
NOC model stems from two factors: (1) the omission of the 
energy derivative term in equation ([8j or the equivalent term 
in equation (|6]i; and (2) the fact that equation ( fTSl i is mani- 
festly non-conservative for neutrino, and consequently lepton, 
number when integrated over mass. In the neutrino opaque re- 
gions, the energy-derivative term is responsible for promoting 
neutrinos in energy as they are compressed, as expected for 
a relativistic Fermi gas and first noted by .Castor (.197Z) and 
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Figure 8. Sequence of electron (solid lines) and lepton (dashed Unes) frac- 
tion profiles for the GR-FullOp (upper panel) and N-ReducOp-NOC (lower 
panel) models showing the formation of the shock and the shock progress 
through bounce into the breakout phase. The colors represent different 
epochs, which are equally spaced in computational time step with time ad- 
vancing from dark to light gray (violet to red in the electronic edition). 

lArnetj ([19771) . 

5. CONTEMPORARY MULTIDIMENSIONAL 
SUPERNOVA MODELING 

5.1. Multidimensional Supernova Codes 

There are five extant codes that can compute the spec- 
tral neutrino radiation hydrodynamics for core-collapse 
supernova simulations in 2D or 3D. These codes are (in 
alphabetical order by name): the FAU-NCSU-Oak Ridge 
2D/3D code CHIMERA (S. W. Bruenn et al., in preparation), 
the Stony Brook 2D code V2D (ISwestv & Mv ra 2009, 2005X 
the MPA-Garching 2D code VERTEX (RamoD & JankJ 
I2002t Burasetal, 2006), the Arizona-Caltec h-Hebrew 
University-Princeton 2D code Vulcan/2D ( Livne et alJ 
l2004t iBurrows et all lIOOl. and the Kyoto-Tokyo 2D/3D 
ZeusH-IDSA code (|SuwaetalJ[2010). 

Of the multi-D codes, CHIMERA and (Prometheus- 
)Vertex include a spherically-symmetric, post-Newtonian 
GR approximation, while the others are strictly Newtonian 
in their grav i tation , hydrodynamics, and neutrino transport. 
iMiiller et al.1 (120101) have updated (CoCoNuT-)Vertex to 
include general relativity in the transport and hydrodynamics 
using the conformally flat approximation. 

CHIMERA, V2D, and Vulcan/2D ti-ansport neutrinos by 
the flux-limited diffusion method (ELD). Vulcan/2D also 
has a non-moment, multi-angle (Sat) mode. Vertex uses 
the variable Eddington tensor (VET) method with a clo- 
sure computed using a spherically averaged, model Boltz- 
mann equation. The ZeusH-IDSA code uses the Isotropic 
Diffusion Source Approximation (IDSA; Liebendorfer et akl 
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Figure 9. Illustration of the "ray-by-ray" transport approximation. The circle 
represents the neutrinosphere and the solid lines represent two independent 
"rays" in the RbR approximation. The dashed lines are tangents to the neu- 
trinosphere and indicate the regions that contribute to the neutrino field at 
points 1 and 2. The "blob" on the neutrinosphere below point 1 is a "hot 
spot" where the temperature is higher than the rest of the neutrinosphere. For 
point 1 , the RbR method will compute the neutrino field as if the entire neu- 
trinosphere has the properties of the hot spot, overestimating the neutrino flux 
and heating. For point 2, the RbR misses the contribution of the hot spot by 
assuming the neutrinosphere properties are only those of the cooler region 
directly below it, underestimating the neutrino flux and heating. 

12009^. which divides the neutrinos into "trapped' and "free- 
streaming" neutrinos, with a diffusion source to connect them. 

Of these codes, only V2D is capable of solving the full 
space-neutrino energy-species coupling of the neutrino trans- 
port that the core-collapse supernova problem requires, while 
all other codes break at least one aspect of that coupling to 
reduce computational costs and simplify code development. 
CHIMERA, Vertex, and Zeus-nlDSA break the non-radial 
(lateral, or angular) spatial coupling through the "ray-by-ray" 
(RbR) approximation, and Vulcan/2D breaks the coupling be- 
tween energy groups and neutrino species. 

In the RbR approximation, the neutrino transport is com- 
puted as a number of independent, spherically symmetric 
problems, referred to as "rays," which allows for the reuse 
of existing ID neutrino transport codes. (See Figure |9] for a 
schematic illustration of the RbR approximation.) RbR meth- 
ods exhibit good parallel scaling for large numbers of inde- 
pendent radial rays, which can be evolved without commu- 
nication while computing the neutrino transport. Typically, 
in RbR codes, the neutrinos in opaque regions are advected 
laterally with the fluid motions and contribute to the pres- 
sure. The independence of the rays artificially sharpens the 
lateral variation in the neutrino luminosity and heating above 
the proto-NS, which results in some regions of the hot man- 
tle being overheated and others underheated. The transport 
studies of Ott et al. (2008) using Vulcan/2D in multi-angle 
mode showed that full multi-D FLD underestimates the lateral 
variation in the neutrino radiation field, whereas RbR codes 
are exp ected to overestimate the lateral variation. iBuras et aTl 
(|2006) concluded from analysis of their RbR models that the 
transient lateral variations in neutrino flux and heating were 
not very likely to have dynamical consequences for the evo- 
lution of their models. The impact of the RbR approxima- 
tion on the simulation outcomes is not precisely known, and 
proper testing will have to wait until one of the RbR codes is 
upgraded to include full lateral transport, as no extant code is 
currently capable of computing in RbR and non-RbR modes 
and there are significant differences between extant RbR and 
non-RbR codes in other respects. 

The authors of Vulcan/2D have chosen to break the energy 
and species coupling rather than the lateral spatial coupling. 
Vulcan/2D implements computational parallelism by solving 
for 2D-spatially-coupled neutrino transport for each energy- 



species group independently, with communication only after 
transport to integrate neutrino heating/cooling from all energy 
groups. The consequence of this design choice is that Vul- 
can/2D cannot easily include either NlS-driven coupling of 
energy groups or the coupling of energy groups through ob- 
server corrections, nor can it utilize more parallel processing 
elements than it has energy-species groups. 

5.2. Opacity Approximations 

CHIMERA and VERTEX include all of the FuUOp opacities 
plus additional corrections for weak magnetism and ion-ion 
correlations. VERTE X also includes t he neutrino-pai r flavor- 
conversion process ( Buras et al.ll2003h . V2D uses the ' Brueriiil 
(1985) opacities, which are similar to ReducOp, but do in- 
clude the energy down-scattering from NES. Vulcan/2D omits 
all of the NIS scatterings in favor of their IS counterparts, as 
does the Zeus-i-lDSA code because energy-coupled scattering 
has not yet been developed for the IDSA transport method. 
Vulcan/2D, V2D, and ZeusH-lDSA use an IPA for EC on nu- 
clei, which cuts off electron capture by nuclei when the mean 
neutron number N > 40, and overestimates it above the cut- 
off, while CHIMERA and VERTEX use the more accurate 
LMSH EC table. 

Some multi-D supernova codes (VERTEX, Vulcan/2D) use 
a single species, i^x = {j^/jt, J^/jt}, to represent all of the 
heavy-lepton flavor neutrinos, while the Zeus-i-lDSA code 
omits them completely. 

5.3. Observer Corrections 

CHIMERA, V2D, and VERTEX include the observer cor- 
rections in the neutrino transport. In the Zeus-i-IDSA code, 
adiabatic compression is properly handled for the trapped 
neutrinos, and 0{v/c) observer corrections are included for 
free-streaming neutrinos. These codes use neutrino trans- 
port based on equation (|3]l, its equivalent to 0{v/c-), or its 
GR equivalent. Only Vulcan/2D neglects the observer correc- 
tions entirely, by solving the neutrino trans port based on equa- 
tion (|7]i. [The transport equation quoted in lLivne et al.l (l2004h 
also omits the ^qv df /dt-teim, which is typically considered 
of ©(v^/c^) and dropped from most 0{v/c) transport solu- 
tions.] 

6. CONCLUSIONS 

We have examined the consequences of removing (1) GR 
effects, (2) non-isoenergetic scattering and detailed nuclear 
EC opacities, and (3) observer corrections from spherically 
symmetric models of core-collapse supernovae. We have 
found that all of these changes, individually and especially 
when taken together, affect the progress of stellar collapse 
and the post-bounce evolution of the shock and core ther- 
modynamic properties in significant ways, in constrast to 
the assessment s made by iBurrows et alJ (120061 12002) and 
iNordhaus et al.l (12010). We have computed variations in the 
shock radius, neutrino luminosities, and neutrino RMS ener- 
gies as large as 60 km, 35 Bethe s^^, and 10 MeV, respec- 
tively, across the four models considered here. 

Omission of GR results in a less compact core and 
an unrealistically more favorable shock progression after 
bounce. Eliminating non-isoenergetic scatterings and simpli- 
fying electron capture on nuclei drastically reduces the core 
deleptonization and expands the homologous core at bounce. 
Omission of the observer corrections dramatically alters core 
deleptionization, the shock position, and neutrino luminosi- 
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ties after bounce, in part resulting from a complete breakdown 
of lepton number conservation. 

The lepton non-conservation and non-promotion of neu- 
trino energy resulting from omitting observer corrections in 
our N-ReducOp-NOC model results in a compact, low-Yj^ 
core and a shock trajectory that is the least favorable of our 
models. The artificial loss of lepton number, lower neu- 
trino luminosities, and the consequent lower neutrino heat- 
ing rate and shallower shock trajectory may explain the lack 
of neutrino-driven explosions in models computed with Vul- 
can/2D (see Burrows et al. 2007), in contrast to the result^jp- 
ported by other s (Marek & Janka 2009; fBruenn et ailBodOi: 
ISuwa et aLllIoTOu Takiwaki et al. 201 1) . 

Moreover, the changes in and Yl, their gradients, and the 
entropy gradients that we see as we traverse the models shown 
here will change the location and strength of convectively un- 
stable regions in the proto-NS and between the proto-NS and 
the shock. The lepton and entropy gradients in the proto- 
NS drive prompt convection, the entropy gradients between 
the proto-NS and the shock drive neutrino-driven convection, 
and these in turn seed and are seeded by the SASl. That is, 
the changes we have documented in this transport study have 
implications for all of the multidimensional phenomena we 
know to be important in multidimensional supernova models 
once spherical symmetry is broken. 

All of the ingredients (l)-(3) above must be included in 
multidimensional simulations of core-collapse supernovae to 
ensure physical fidelity. Their omission is not the only ap- 
proximation used in current multidimensional simulations, 
some of which (like the ray-by-ray approximation) are inade- 
quately understood and need to be better understood or phased 
out. Certainly, further examination of these approximations is 
required within the context of multidimensional simulations. 
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